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Preface 

The science of statistics can bo of groat use to goodosists and 
photogrammetrists in reducing and analyzing data obtained from meas- 
urements. This thesis will oxploro the applications of statistical 
hypothesis tests in the fields of geodesy. 

Hie data for examplos have boon taken from various studios in 
the field of geodotic science. In each case tho sourco of tho data 
is noted. Often tho conclusions drawn by statistical tests will not 
agree with the conclusions drawn by tho original oxperimontor. It 
is not tho intent of this thesis to criticize tho work of others, 
but only to give examples of how statistical tests may be usod to guide 
geodesists in drawing conclusions from obsorved data. 

The author gratefully acknowledges the guidance of Dr. Urho A. 
Uotila and the assistance and inspiration of his wife, Judy, in the 
preparation of this thesis. 
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CHAPTHt 1 
INTRODUCTION 

In geodesy and photogrammotry, large quantities of data are 
col looted in the form of monpuromontrj such no angles, lengths, and 
gravity values. It is imperative that the relevant information con- 
tained in a mass of geodetic data bo expressed by comparatively few 
values. To accomplish this task, geodesists, as scientists in other 
fields of physical and social science, have found many solutions through 
the science of statistics. 

Statistics is the study of populations, variances, and methods 
of data reduction. Some statistical methods, notably the method of 
least squares, have found universal acceptance in geodosy. Other 
statistical methods are not so widely used. It is the purpose of this 
thesis to exploro tho applications of statistical hypothesis testing to 
the problems of geodesy, 

1,1 Statistical Theory 

A STATISTIC is a value calculated from an observed sample with 
a view to characterizing the population from which it is drawn. 

Statistics which will bo of value to geodesy are: 

Population mean yH 

r» 

— X * 

Sample mean X 3 ^ — - 



True variance 
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z 
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Variance estimate^ 

Fisher's F Statistic 
Pearsen’s *X. 2 Statistic 
Student's t Statistic 
The standard error 
Hie estimate of standard error 1 s 
A statistic which, on the average, gives the right answer is 
said to be unbiased . If a statistic gives values which are concen- 
trated more closely to the right value, the statistic is said to be 
efficient . 

A widely used concept in statistical theory is a random 
variable, A random variable is a quantity which takes on a definite 
value at every point of a sample space. Geodetic measurements are con- 
sidered to be random variables, 

1.2 Density and Distribution Functions 

Assume that a sample space is such that each point of the sample 
space can be characterized by the value of a continuous variable, x, 
which can take on all values between -co and +co . If the event, A is 
defined as the set of all sample points characterized by the inequality 

x ^ x° 

where x° is some fixed value. The cumulative probability distribution 

function, F(x°), and the probability density function, ^(x) , are de- 

iThe letter m is commonly used in geodetic literature. In most statis- 
tical literature Greek letters represent true values and the Roman equiva- 
lent represents the statistical estimate of this value. For ease of 
understanding, the statistical convention will be followed. 

^See discussion by J.L. Steam (1964). 
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fined by 



P(A) - F(x°) 




— CO 



The density function, j2f(x) , has the additional significance that 

P(x°^ x x° 4* dx) « j0Xx°)dx. 

Since the total probability for a sample space must equal unity, 
must include a normalization factor 



F(oo ) 



1 



-f co 

^(x)dx 



1 . 



The basic requirement for any probability density function is 

that this integral exists, A further requirement is that 

jtfix) ^ 0 for - 430 < x < + . 

There are many distributions and density functions. As an 

example the uniform density function is defined 

^(x) - 1_ for - a x -f a 
2a 



and 



jtf(x) - 0 otherwise. 

Figure 1 shows the corresponding probability density function, 
^(x), and the cumulative probability distribution function, F(x) , for 
this uniform density function. 





Figure 1 
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A quantity is said to bo normally distributed when it takes all 
values from -<r> to +<o , with frequencies given by a definite mathe- 
matical law, namoly, tho logarithm of the froquency at any distance, d, 
from tho centor of tho di stribntion is less than the logarithm of the 
froquency at the ccntor by a quantity proportional to^ 2 . Tho distri- 
bution is symmetric with greatest frequency at the center (Fisher, 1925)* 
The density function of a normal distribution of moan m and variance <?* 
is given by tho expression 

, , _ (x 

e . 

Figure 2 shows the probability density function and tho cumu- 
lative probability distribution function for tho normal distribution. 

The scale of x can be changed by measuring each x value by its distance 
from the moan, and adopting tho standard deviation 6 as a unit of meas- 
urement. An ordinate of this normal curve is then 
W * 

6 ~ 

The quantity w is callod the standard normal deviate (Handel, 1964), 

A unit normal deviate is a distribution which has a moan of 0, and a 
variance of unity. 





Figure 2 
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The normal distribution is very useful in statistics because of 
the Central Limit Theorem. This theorem may be expressed as follows: 

Given a population of values with a finite variance, if inde- 
pendent samples are taken from this population, all of size N, then the 
population formed by the averages of these samples will tend to have a 
normal distribution, regardless of what the distribution of the origi- 
nal population is; the larger N, the greater will be this tendency 
towards normality. 

Thus far only functions of one variable have been discussed. 

If there are more than one variable associated with each point in sample 
space, multivariant functions may be defined. For example, a multi- 
variant probability density function is defined as 

^(Xi 0 , x 2 °, x rt 0 )dx 1 dx 2 dx n 

- Ptxj 0 ^. 6 . Xi° + dxj , i - 1, 2,..n). 
The right side of this equation may be interpreted as the probability 
that all the inequalities hold simultaneously. For an excellent discus- 
sion, the reader is referred to pages 16 through 18 of Hamilton (1964). 

1.3 Statistical Hypotheses 

Webster defines a hypothesis as a tentative theory or suppo- 
sition provisionally adopted to explain certain facts and to guide in 
the investigation of others (Webster, 1954). A statistical hypothesis 
is thus a theory about some population. 

The only way that one can be absolutely certain of the truth or 
falsity of a statistical hypothesis is to examine the entire population. 

Since measurements can take on an infinite number of values, exami- 
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nation of tho cntiro population in gcodotic -applications is impossible. 

One is thon forced to mako a decision basod on a fow measurements. 

Statistically speaking those measurements are a sample taken from tho 

population. Tho procoss of using this samplo to test tho truth or 

falsity of a hypothorjin is called statistical tests. Thero is in those 

tests no cortainty that a mistako has not been made. There aro, in fact, 

two different kinds of orrors which can bo made. These aro called ; 

Type I (& ) error — tho rejection of a hypothesis 

which is truo 

Type II error — tho acceptance of a hypothesis 

which is falso 

These errors will be discussod in detail lator. 



Table 1 

Definition of tho Types of Errors Associated 



Decision 


Truo Situation 


Hypothesis is Truo 


Hypothesis is False 


Accopt tho 
Hypothesis 


No Error 


Typo II Error 


Reject tho 
fiypothesis 


TVpe I Error 


No Error 
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CHAFUSR 2 

HYPOTHESIS TESTIS 

It may bo of valuo to know if data are normally distributed. 

A simple test is to compare tho histogram to a normal curve (Dixon, 

1957). Tho percentage of tho data in a given group of tho histogram 
can be compared to the area under tho normal curve corresponding to the 
given group. 

For example; givon the mean of a data sot as 30* and its stand- 
ard deviation of 5* 80$ of tho data is found botwoen 20 and 35* Is 

this data normally distributed? 

The standard normal deviates for tho data group are 

*20 = 20-30 = - 2.0 

x 35 = 35-30 =1.0. 

From Appendix I the area under the normal curve is 

from -CD to 1.0 = .8413 

from -OO to -2 e 1.0 - .9772 = .0228 

The aroa under the normal curvo = .8185 

or 81.85$. This would indicate that the data in this group has a distri- 
bution close to a normal distribution. 

The °)t z statistic, discussed in section 2,4, may also be used to 
test the distribution of a group of observations. A discussion of this 
test will be deferred until the statistic has boon introduced. 

For an example of more elaborate tests for normal distribution 
the reader is referred to the paper by Steam (1964). 
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2.1 Tests Involving the Normal Distribution 

If a single observation, x, is made from a normally distributed 

2 

population of mean,^, and variance, G , statistical theory states 
that (Hamilton, 1964) 

(1) w - (x -M) 

S' 

has a normal distribution. The probability that the magnitude of w 
calculated in this manner exceeds some specified value 

(2) P(l w/ > w r ) 

where is the value of w to the right of which lies an area jr under 
the probability density curve is 

(3) P( |w| > w ff ) - # . 



We can then write 

(4) P(|w/ > w,) - _J_ 



Yz'i 



f 



dw + 



+ CD t -i 



From Figure 3 we can interpret the value of the integral as the cross 




Figure 3 

Mathematically Hamilton (1964), shows 



(5) 



P(/w| > w^) - 2F( - w y ) 

- 2(1 - F(wj.)) 

- 2 V 



The value of F(w^) can be found from a table of the cumulative normal 
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distribution function (Appendix l) # 

If is 1*96 

PCI w/ ;> 1.96) - 2(1 - 0.975) - 0.05. 
Equation (5) can be expressed as 

( 6 ) 



or 



P(| x -M l > » 20- - F(vy)) 



(7) P(^-6-w^< x^^ +<5wp - 1 - 2[l - F(w^ 

- 2F(w r ) - 1. 



For the example 

P (/(- 1 . 96 <s x^/{ + 1 . 96<r) - 0.95 

that is, the probability that a single observation, in the normal popu- 
lation given, will lie within 1.966“ of the mean, is 957** 

The inequality expressed in (6) forms the confidence interval 
which is the basis for the test of a hypothesis involving the mean of a 
normal population. 

Symbolically this null hypothesis, H 0 , can be expressed 
(8) H 0 : A -/(o 

and the alternative 



Hi 'A\Ao- 

The unit normal deviate is calculated from (1), if | w| >1.96 
the hypothesis, ll 0 , can be rejected. The risk of rejecting a true H Q 
(TYPE I ERROR), or the level of significance Cck ) is .05. It should be 
noted that any other level of significance could be selected simply by 
selecting a different value for from the table. As an example, if 
w were 1. 64, c< would be .10. 

The areay, shown in Figure 3, varies directly aso( . The 
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shadod area, 2 6 is in this case equal to c* . This region is known as 
the critical rorri-on (Ostle, 1963). l f the value of tho tost statistic 
used in a particular test of a statistical hypothosis falls in this 
region, the hypothosis is rejocted. 

If oc= 2 y , the test is said to be two tailed : that is, the 
hypothesis was rejected either if 

w > 1.96 
or 

w< -I.96, 

The test can also be formulated such that the hypothesis is re- 
jected only if 

w > 1,96 

or if 

w< -1,96, 

In each of those cases tho probability of a Type I orror (rejection of 
a true hypothesis) is tho area under only one tail of the curve and 

( 9 ) = Y 

Figure ( 4 ) shows the critical region for the one-tailed test. 




Figure 4 
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A tost of a statistical hypothosis based on either extreme of 
tho distribution is a one tailed test . It is interesting to note that 
for ar\y symmetrical distribution a one tailed tost at 100 (cx ) $ signif- 
icance level is equivalent to a two tailed tost at 100 ( 2 <*)^ signifi- 
cance lovol. This is true since the same critical value of w applies 
to the two cases. 

2,11 Test Procedure 

In order to test the hypothesis 

H 0 : *=M 0 0 = 

against the alternative 

H 1 : *=* 0 

at a significance leveled, tho following steps would bo followed: 

(1) Select a lovel of significance oc, and find the corre- 
sponding Wy from tho appropriate table. 

(2) Compute tho unit normal deviate 

w = * 

CT 

( 3 ) Re.iect the hypothesis if lw|>wy . 

For example, a set of observations has a sample mean X of 1.50* and a 

variance, of 0.25. Could tho true mean,// , be 2.00? The hypothe- 

^ o 

sis to be tested is 

H Q 2.00 

If oc is selected to be 0.05* the tabulated value of w g is I. 96 . 

w = 1.50 - 2.00 = ^ = -1 

.5 .5 

1 < 1 .96 

Therefore, the hypothesis can be accepted at the level of significance 
selected. 
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The testing of a hypothesis such as this has very little 

2 

practical significance in geodesy, since the variance, <5 , is seldom 
known. 

This test has been discussed in great detail because the prin- 
ciples involved apply to all statistical tests, regardless of the 
statistic, or distribution used. 

2.2 Testing with a Sample Mean 

The sample mean of a set of observations is defined as 

X - I H x i 

n i-1 

where n is the number of observations. 

Statistical theory states that the sum of normally distributed 
random variables is normally distributed (Hamilton, 1964). 

The density function of X is 

& (X) - 1 f_n ^exp -n(X -^y) 2 . 

6 (2 rrj 2S 2 

The standardized variable is then 

w - (X -x#)n^ 

< 5 " 

Following the same method used for an individual observation, the proba- 
bility that the deviation of the sample mean exceeds a specific value is 

P(/(- 6~y r zX^yt(+ ^ r ) - 2F(w / ) -1 
n 2 n2 

To test a hypothesis from a sample of n observations and a sample mean 
of X, 

H 0 : H-// q 

against 

Hi -M\M Q 

compute 

, w - (Y - 



er 
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The hypothesis is rejected at the 100c*7<» level if |w/ > w^. 

With this test a cpnfidence interval is established around X. 
The probability that the true mean,^, lies within this interval is 



vations is part of an established population with a mean/(, and a 
standard error of 6^. 



the level of significance chosen for the test. There is also the 
possibility of committing a Type II error, that is the error of accept- 
ing a false hypothesis. 

The probability,^?, of a Type II error is dependent upon the 
specific alternative hypothesis which is presumed to be true 



when is true* From (1) 

w • (X 

_ 

If the alternate hypothesis is true, w^ is distributed as the 
unit normal deviate, and 




This test would be useful to determine if a new set of obser- 



It has been shown that the probability of a Type I error is 



Hi xA-H, 



The probability of a Type II error is 



( 10 ) 






(ID 



W 1 “ n^(X -Ap *Aq 




This can be written as 



( 12 ) 




Equation (10) can then be written as 



6 ~ 

written as 



( 13 ) 




0 ~ 
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(14) /? - P( 




It should be noted that as ) becomes smaller the proba- 



bility of not rejecting a hypothesis when it is false is nearly as great 
as not rejecting it if it is true. The power of a statistical test is 
defined ad 1 The power of a test increases as ~/J t becomes 



larger or, as the number of observations increases (Hamilton, 1964). 
Figure 5 shows the power of the test for c* .05 and .01, as a function 



Figure 5 

Note that 1-/5 varies directly as The value of is 
determined by the value chosen for c^, If the critical region is in- 
creased, decreases. Concerning this problem Graybill (1961) states. 



"We should like to minimize the likelyhood, or proba- 
bility, of making either of these two errors. However, 
in general, for a fixed number of observations, if we 
decrease the probability of making an error of one type, 
we increase the chances of making the other." 

Power curves can be found in Table A-12 of Dixon (1957). 



of c f where <f - (#1 -/Qn^ . 




-S.o 



o 



+ S.o 
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Using the t distribution the tost 

against 

Hi •■A*A 0 

can be carried out. 

First compute 

t = (X - M q W = (X -/( Q ) . 

s o S X 

Reject the hypothesis at the cx, levol of significance if 

| tl > *71-1 ,<* . 

This test is valuable to the geodesist attempting to deter- 
mine if a new set of observations is from the same population as 

— 2 
previous observations with a mean of X and a variance estimate s . 

This problem often arises in the weighting of observations in an adjust- 
ment. 

The following problem will serve as an example. 

In his evaluation of the Laser Theodolite, Dunn (1966) 
observed the horizontal angles botwoen two targets at various current 
levels. At 35 amps the angles obtained were (Dunn, p 52): 
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3° 


Loft 

17' 22". 2 


3° 


Right 

06 * 29". 8 


3° 


Moan 

06 • 56". 0 


Residual 

-5.0 




07 


17 .5 




06 


36 .8 




06 


57 .1 


-3.9 




07 


00 .3 




06 


43 .6 




06 


52 .0 


-9.0 




08 


29 .2 




06 


39 .8 




07 


34 .5 


33.5 




07 


18 .4 




06 


33 .8 




06 


56 .1 


-4.9 




07 


13 .1 




06 


27 .8 




06 


50 .4 


-10.6 


mean 
3° 07' 


26". 8 


3° 


06 ' 


35". 3 


3° 


07' 


01". 0 




l™3 


= 1379.83 


£vv J _ 

n -f 


s 2 = 276.0 




5 0 


= + 16" 


.5 



lh© t tost may bo used to compare tho moan of ono set with a 
grand moan of previous sots. 

From th© previous two sots of observations th© moan is 
20 M amps 3° 06' 51". 1 

30 M amps 3° 06* 52". 1 

mean of 1 & 2 3° 06 ’ 51" .6 

Tost th© hypothesis 

H 0 X M 35 = 30 06 * 51 ". 6 

against 

Hi : M. 35 £ 3° 06* 51". 6 

t = (3° 07' 01". 0 - 3° 06 ' 51 ". 6 ) ( V~6) 

t = 9 ".4 (2.45) = 0.990 
16.5 

^5» 5$ = 2 *57. 

Sine© t < tj this hypothesis can b© accoptod at th© 5$ signifi- 



cance level 
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2,4 Tests Involving Variance 

The statistic ^ was introduced by Pearson. It is defined by tho 
density function 



(17) 



- *>2 



2^T(%) 



{x z f /z 1 



for ~K 3 > O othorwiso <f) ( "X 2 ) = O . 

It is also known that the density function for tho variance estimate 
s c from a sample sizo n of a population with a true variance 6 is 



(Hamilton, 1964) 
for s 3 > O 



(■ 









T(y/2.) I 2 6'V < ^ X P['Y§ t -] Cs ^ 

otherwise <$ ( s* ) =0. 



v/z 



Setting X*- these two donsity functions are identical, Ihe 

value of is distributed as 
<s 

The value of X^uch that 

= « 

are tabulated in Anpendix III. 

The*XL 2 statistic can be used to test whother the true variance 
estimated by s^ is equal to some variance G 0 Z . Tho hypothesis can be 
stated 

H 0 * e 2 = <T e 2 
: G^ \ Gj , 

To test, compute C X. 2 = 

°o 

The hypothesis is rejected if 



“X. < °x‘ 






or 






Y>°Vz 
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The uso of this statistic to test hypothosis concerning 
geodetic and photogrnmmotric data is vory limitod. To computo 
tho statistic, G x tho truo variance of the population must be known. 

This quantity is rarely known for tho typos of data analyzed in geodesy. 

One application of the Chi-squaro tost would be in trian- 
gulation. Tho square of tho dosirod standard error of the not could 
be considered as tho truo variance, S'*. Tho Chi-squaro statistic com- 
puted at each station in tho net, at the time of observation, could 
then be used as a criterion for acceptance or rejection of the obser- 
vations. 

The "X 2 statistic can bo used as a tost of distribution. To 
test tho distribution, observations are grouped into n groups accord- 
ing to their values. The expected number of observations in each 
group is computed on tho basis of the assumed distribution. Let n^ be 
tho number of observations actually found in tho kth group and d^ bo 
the number predicted by the assumed distribution, Tho statistic 

‘X 2 * Z -‘M* 

k’l d* 

is distributed approximately asl^_„if each d^ is at least 5 (Hamilton, 
1964). 

2.5 The F Statistic 

Tho F statistic or varianco ratio, first introduced by Fisher, 
is a useful statistic. 

F is defined 

(18) Fv 1^2 " (Y 1^1> /(Y9Ayo) 

where Y^ is distributed as *X ? with \p dogroos of freodom, and Y2,‘ 

INDEPENDENT of Yj, is distributed as ^with ^ dogreos of freedom. 
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The probability density function is given as 

0(f) * T(rv?, ^pj/2} (& Y ¥f)~ ‘ Wt 

~r (y/ 2 ) T(vji) \ v, / r V v z / 



for 



F 0 



othorwi3o 



a 

0(F) =0. 



It was shown previously that 

2 



X 



V5\ S." 



For two samples from two normal populations 



nr 

<r,‘ F ' 



3 1 

So 2 



is distributed as 



Ihus the hypothesis that tho samplos wero drawn from populations with 



identical variancos 



H 



o 




=<r 2 



2 



can be tested against 

% s 6 i 2 ^ 2 2 

compute 

F = s 1 /s 2 . 

The null hypothesis is rejected if the computed statistic falls 
in the seloctod critical rogion at either ond of the distribution 

F > F ^l, \? 2, 1 -o^/2 or 

F ^ \? 2 , oc/2 (Dixon, 1957). 

The inequalities shown above dopond upon tho definition of oC , 
In some literature they will appear reversod. 

Tables of the percontiles of tho F distribution can be found 
in handbooks and statistics texts such as Dixon and Massey, Table A-7c. 
Appendix IV shows a sample of such a table. 

As an example of the F test, the data takon by Dunn (1966) and 



Sprineky can be analyzed 
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Lt, Dunn made 3^ observations of a horizontal angle with tho 



Laser Thoodolito. 



The values were 



3 06' 56". 1 
06 44 .7 
06 50 .4 
06 51 .6 
06 47 .2 

06 50 .6 
06 56 .2 
06 4? .0 
06 53 .9 
06 51 .1 
06 54 .6 
06 47 .7 
06 40 .5 
06 37 .7 
06 40 .2 
06 45 .5 
06 45 .8 
06 49 .7 
06 56 .0 
06 57 .1 
06 52 .0 
06 56 .1 
06 50 .4 
06 35 .0 

06 22 .4 

07 00 .1 

06 58 .4 

07 04 .3 

06 45 .3 
06 19 .6 

07 00 .7 

06 34 .8 

07 17 .6 
06 18 .2 



Mean of 34 observations 3° 06* 48". 3 
Standard deviation + 12", 4 

153.8 



Sample variance 
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Using the same equipment Capt. Sprinsky obtained the following 

5° 27' 19". 9 

26 58 .4 

27 00 .0 

26 53 .4 

26 54 .3 

26 57 .1 

27 13 .7 
27 22 .6 
26 58 .2 

26 41 .0 

27 07 . 1 

27 13 .2 
26 48 . 1 
26 41 .4 
26 30 .7 
26 44- .5 
26 49 .5 
26 45 .0 
26 52 .2 
26 33 .2 
26 38 .8 
26 39 .2 

26 38 .0 
26 28 .1 
26 40 .2 

26 29 .1 

27 10 .7 
26 42 .7 

Mean of 28 observations 5° 26' 50" .37 
Standard deviation + 14" . 3 

Sample variance 204.5 

Were the true variances of those two sots of observations equal? 
As a hypothesis this question can be stated 

H o :< ^l “6~2 



and can be tested against the alternative 

Hi : *Y < d 2 



The test statistic is computed 



F = 204.5 = 1.33 
153.8 
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From a tablo of porcontilos of F distribution, such as tho sample in 
Appendix IV, the values of F aro found to bo 

f 28, 35. 2.5# = * 469 
f 28, 35. 97.5# = 2 ‘ 01 ' 

Since 

F >f 28, 35. 2.5# 

and 

F 4 F 28, 35. 97.5# , 

\*e can accopt the hypothesis at tho 5# significance lovol. It can then 
be concludod that tho observations of Lt. Dunn and Capt. Sprinsky have 
tho same true varianco. 

The following example will demonstrate another uso of the F 

test. 

In his study of tho Kern DKM 3 thoodolite Abby (1965) moasurod 



angles with a Wild T-3» fl nd tho DKM 3 using the center wire and the 5 
wire field. Thoso rosults wore obtainod: 





Instrument 


Observations 


Wires 


s 


s 2 


1 


T-3 


16 


1 


0".99 


0.980 


2 


T-3 


16 


1 


.78 


.608 


3 


DKM 3 


15 


5 




.221 


4 


DKM 3 


10 


5 


.44 


.194 


5 


DKM 3 


10 


1 


.66 


.436 


6 


DKM 3 


30 


5 


.61 


.372 




Tho F 


statistic can bo used 


to assist in drawing 


conclusions 



(a) Are the 5 wire observations significantly hotter than the 
1 wire observations with tho DKM 3? Tho hypothosis to be tested is 
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H 

For this examplo a one-tailed 
can be rejected if 



: <5 






o * ' J L 
tost will bo used. 



The null hypothesis 



sj 2 \/Tl -<* ) (ni - 1) (n 2 - 1). 

i? 

Comparing sots 3 and 5 



F = 0.4-3 6 = 1.9 73 

0.221 

F (l^, 9) 5$ = 3.00 

The difference is not significant at the 5$ level. 

Using sets 4 and 5 

F = 0,426 = 2.247 

0H& 

F (9 . 9) = 3.18 

The difference is not significant at the 5$ level. 

And using sets 5 and 6 

F = 0.436 = 1.172 
0.372 

F (15. 29) 5# = 2.03. 

Again the difference is not significant at the 5$ level. 

(b) Are the 5 wire, 30 observation sets significantly better 
than the 5 wire 10 observation sots? 

F = 0.372 = 1.918 

0.I9E 

F (29, 9) 5 $ = 2.86 
Ihe difference is not significant. 

(c) Is the DKM 3 significantly better than set 1 with the 



T-3? 
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F = ,980 = 2.24? 

'.436 

f (15. 15) % = 2M 

Again tho test shows that the difforonco is not significant, 

(d) Is the DKM 3 . 5 wire procedure significantly better than 

the T-3? 

F = 0,980 = 4.46 and F = ,608 = 2.74 

0.221 .221 

F (15. 14) = 2 * 44 

The tost shows that the DKM 3» 5 wire procodure is significantly better 
than the T-3. 

This problem also shows that statistical tests are not in- 
fallible. From parts (a) and (b) one can conclude that the dif- 
ference in procedure is not significant. Part (c) draws the con- 
clusion that the 1 wire DKM 3 procedure and the T-3 are not sig- 
nificantly different, yet (d) concludes that the 5 wire procedure 
is significantly different. 

Statistical inference must bo used with caution. Judgement 
must also be used in drawing conclusions. 

2.6 Testing Correlations 

The distribution of the sample correlation coefficient, r, 
is very complicated. Often in geodesy it is valuable to know if 
the true correlation coefficient, p , is zero. If p is equal to 




zero the statistic 



2.6 



is distributed as Studont's t with (n-2) degroos of freedom (Guttman, 

1965 ). 

To test the hypothesis 

H 0 ;/ = 0 

against 

Hi : j> ± 0 

compute the statistic 

C = r Vn-2 

TUrZ 

If C > tn-2, e</2, hypothesis may be rejected in favor of the 
alternate. 

To test the hypothesis 

H 0 : f> = f, , where 

we must know the distribution of P when -j 3 £ 0. In most geodetic 
applications it is more important to know if correlation exists, 
rather than if it might be some specific value. For a discussion 
of the distribution of P £ 0 the reader is referred to Graybill (I 96 I). 
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CHAPTER 3 

TESTS INVOLVING MULTIPLE VARIABLES 

So far the tests which havo been discussed doal only with one 
variablo. The techniques of Rogrossion Analysis allow the testing 
of hypothesis involving two or more related variables. 

3.1 Regression Analysis 

The functional relationship between some obsorvable L and 
variables can be expressed mathematically 

(19) L = f(Xi, X 2 X n /€?, e.) 

g> £ - is a parameter in the function. Hiis equation is ofton abbre- 
viated 

L = f(X if X 2 X n ). 

To the statistician this equation is known as a regression 
function. Ihe geodosist and photogrammetrist will recognize it as 
the mathematical structure for a method of observation equations 
adjustment problem. 

Tho mathematical structure of tho problem may bo chosen 
by two methods. In goodosy the analytical consideration of the 
phenomenon involved is tho preferred method. Tho examination of 
scatter diagrams plotted from the observod data can also yield a 
workable structure. 

3.2 Fitting a Linear Mathematical Structure by Least Squares 

Suppose that a linear relationship exists between a dependent 
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variablo, y, and an independent variable, x, such a3 tho relation- 
ship between gravity and height in froe air. From observations n 

sots of points, (xj,yj), (x 2 , y 2 ) (x^, y n ) aro obtained. Tho 

sample structure can bo expressed 



y = /o +f 3 \ x 

From tho obsorvod data wo wish to obtain the values of the unknown 
/?o,/?i and 

Let us assume that y for any given value of x is a random variablo 
which is normally distributed with a mean x and variance 6~ 2 . 

We will also assume that /?, and /3, do not depend upon x , The con- 
ditional probability density function of y for a given x can be written 

(20) f(y/x) = 1 exp f -1 (y -/g '- A x) 2 7. 

Wari- L E<s - 2 J 

The expected value of y, E(y/x) is 

(21) E(y/x) =/ D + /^x. 

The conditional probability density function^ of a set of n obser- 
vations (x,y), (x 2 , y 2 ) (xjj, y ) is 

(22) TT-f (V%.) = jr exr[-££i ‘A *0] 



{&*'f -a -a * fjj 



Applying the principles of maximum likelyhood (Guttman, 1965) the 
parameters which maximize tho probability density function must be 
found. To maximize this function tho quantity 

(23) — 2 



s - £ (Yi -/^ 0 "AV 



3 



For proof see Guttman (1965), Appendix III 
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must be minimized with respect to and ^5, . The process of esti- 

mating !& 0 and by minimizing the sum of the squares of the residuals 
is tho well known method of least squares . Tho value of and ^3^ 
which minimize F (/^,/f) aro thoso for which 



(24) 



ZF(&A) 




=0 

= 0 . 



It should bo notod that for a normal distribution, tho method of 
least squares gives a maximum likelyhood estimato of the paramotors. 

If the error in tho observed quantities is not normally distributed, 
the method of least squares will give a minimum variance estimate of the 
parameters but this estimate will not bo tho maximum likelyhood estimate. 
Taking the derivitivos of (23) we obtain 
(25) -it (y t vsyp - 0 

-2^x i (y i a -/?iXj) - 0 

c-i 

These equations can be rewritten as 

T % x i 

x iyt £ x i +/^i ^ x i 2 

In this form they aro referred to as tho normal equations. In matrix 
form this equation is the formula 

NX + U = 0. 

If the determinant of N is not zero a solution to this equation exists. 
The values are 



b o 




b L x 



b l “ (x t - x)(y i - 7) 

(xi - ~k) z 

whore b Q and bj are estimates of jQ o and jg (Guttman, 1965). 
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Graphically all valuos of x and y could bo represented as a 
straight line in the x, y piano with intercept b Q , and slope b^. This 
lino is called tho regression line. 

The unknown parameters can bo found using tho matrix oquation 

NX = -U 

or X = -N-^U 

tho variance of unit weight, s G 2 , (or m 0 ^, as frequently usod in 

geodetic literature) is given by 

(28) s 0 2 = V'PV 

n-u 



In tho linear regression problem u is 2. From the matrix adjustment 
we also form the weight coefficient matrix, Q x . Tho variance - 
covariance matrix Z, is formed by multiplying Q x by s 0 2 . In tho linear 
regression 



(29) 



’b o 



S bb. S b 



z 



Those quantities are the statistics which will be used in tho testing 
of hypotheses concerning this regression. 



3,3 Tests of Hypotheses in a Linear Regression 

In most linear regressions the estimator of greatest impor- 
tance is the slope of tho line, b^. To tost if the slope is signifi- 
cantly different from some hypothesized value, say^ the t statistic is 
used. Stating the hypothesis 



H . : A-A' 



the alternate is 
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the statistic computed is 

(3 0) t = (b t -/? 3 ') 

s . 

The hypothesis is rejected if 

1 * *(1- / 2) (n-2) 

or 

t -" t (l-/2) (n-2) 

It should be noted that if y is independent of x, will be 0. We 
can test for this by setting ^ oqual to 0, then 

(31) t; = b- - 0 

S iV 

would be the test statistic. 

Other hypotheses concerning the linear regression are 
O) l <0 =/*<>= <3 o’ 

(2) H o 'A =A' and 

For hypothesis (l) the t statistic is used 

(32) t = (b 0 

Sb o 

The F statistic can be used to simultaneously test hypothesis 
(2). This method will be discussed in soction 3.31. 
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Table 2 

Summary of Tost Procedures 



Hypothosis 


Statistic 


Equation 


Rejection Region 


A. 'A' 


t 


30 


|tl 2^ 4 t ,-«v a j 


a.' a: 


t 


32 


|t| > 


A -A' 

and 


F 


34 

(soc. 3.31) 





Tho mathematical structure discussed is quite general. For 

example 



y = A* + A sin t 

can bo handled by the methods previously discussod simply by sub- 
stituting 

x = sin t. 

The expotontial problem u = 16 o^ v reduces to a linear form 

log u = logV'scTv. 

Making the following substitutions 

y = log u 
£, = log y 

X = V 

Wo have tho problem in tho form 

y -/? 0 +^ L x (Ostle, 1963 ). 

As one example of the goodetic applications of a linear 
regression data presented by Laurila (1965) can bo analyzed. In this 
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report , Dr. Laurila assumes that tho relationship between absolute 
humidity and measuring error is linear. The mathematical structure used 
is 

D t - D 0 = CA + K 

where D is tho observed distanco, D+ the true distance of the reference 

o u 

base, A denotes tho absolute humidity, C represents tho humidity 
coefficient, and K is an instrument constant. 

From a least squares adjustment of 25 measurements the following 
results were obtained (Laurila, 1965)# 



c 


s c 


K 


S k 


0.97 


+0.10 


-11.4" 


+1.2 



Figure 6 is a scatter diagram plotted from tho data given in the 
report, for tho reverse readings. 

Using tho t test, hypotheses about this lino can be tested. It 
should be kept in mind that these tests are conditional tests. They 
take into account only the effect of errors of the parameter tested. 



3^ 




V£> 



Vf> 



fa) 






\o 

i 



O 



R 8 SOL UT£ HUMIDITY C 3 r / m 3 ) 



( 1 ) 
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H 0 : C = 0 

t = 0.22 = 9.7 

.10 

*24, ~ 2 * 0 ^ 

This hypothosis can bo rejected. 

(2) H 0 : K = 0 

t = 11.4 = 9.5 

1.2 

*24, 5 ^ = 2 *° 6 

This hypothesis can also be rejected. 

From a combined adjustment (Laurila, 1965 ) the parameters wore 
found to be 

C = 0.66 
K = -6.8. 

Were these the true parameters for the lino given by the forward readings 
only? 

(3) H : C = 0.66 

o 

t = .97 - .66 = .31 = 3.1 

.10 .10 

t 24, 2 *° 6 

This hypothesis may also bo rojocted, 

(4) H q : K = -6.8 

t = 11 .4 - 6.8 = 4.6 = 3.84 

1.2 1.2 

Again this hypothesis may bo rojocted. 

In a similar manner the linear coefficients obtained from a 



least squares adjustment may bo tested for a zero value or some 
theoretical or previously determined value. 
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3,31 Matrix Approach 

The simplest way to handle adjustment computations is to use 
matrix algebra. It would be valuable to have matrix methods to test 
statistical hypothesis concerning the results of an adjustment. Hamilton 
(1964) has given such methods. 

From an adjustment computation, using the methods outlined by 
Dr. Uotila fiv66) the variance-covariance matrix^!, the solution matrix 
X and the estimated standard error of one observation of unit weight 
s Q ^ can be obtained. From the derivation 

X x - So 2 QX - s Q 2 ir l 
Hamilton (1964) shows that 

(X - XJ’Z " 1 (X - X) 

is distributed as 

^ u^ u 3 ^u, n - u 

*X 2 n - u/(« > u> 



by matrix algebra 



Applying Hamilton’s test 



s 2 Q 




- 1 



(N). 

~T~ 



H c : X - X H 

where X^ is the hypothesized value of X, compute 
(34) S H - 1 (X - XjP' Q* 1 (X -%) 

“ " T 



If Sjj exceeds F Uj n . u^^he hypothesis, H 0 , may be rejected at the o< 



XT 



level of significance. 

Hamilton (1964) also developes the theory needed to test 
hypothesis when constraints have been placed on the structure. 
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Introducing Lagrange multipliers, K / ,the function 

0 = V'PV - 2K' (CX + Z) 
is obtained. Differentiating this expression 

d0 = 2V'PdV - 2K'CdX. 

Let X denote the least squares estimate under the constraint 

d0 = 2 £-U + NX - K'Cj dX. 

Minimizing the residuals 

0 = 2 [-U + X'N - K'C JdX 
so 



Substituting 



K'C = -U' + X'N 



O' = X*'N 



where X* is the best least squares solution without conditions,, the equa- 
tion becomes 

K'C = (X - X*)'N 

KCN” 1 c' = xc * “ X*' C' = Z' - X*C' 

thus 

K =(Z' - X*C' ) (C*N -1 C' ) _1 

eliminating K from the preceding two equations we obtain 

(X - X*)'N = (Z - X*C') (CN“ 1 C , )“ 1 C 
or 

(35) X' = X*' + (Z - X*'C») (cn” 1 c j^cn" 1 

The weighted sun of the squares of the residuals is given by 

Rq = V’PV + (X - X*) , (A f PA) (X - X*) 

V is the matrix of residuals without the constraints. 
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The expected value of Rq ic (n - u + b) whore b is the numbor of 
conditions. The expression may be rewritten as 

= ” R o 

whore 

R c = V'P V 

That is, R Q is tho unconditional least squares sum. and R^ is the addi- 
tional sum of squares due to the constraints. Hamilton (1964) shows 
that 

R H = (Z - CX+)» (C(A'PA)" 1 CT 1 (Z - CX*). 

The ratio 



is distributed as 



% = r Q - R o 
R o R o 



n - u 



Substituting tho values of R^ and R 



-w-1 



Ik = (z - cx*)» (c(a i pa)~ 1 c i ) (z - cx*) . 



R, 



(n - u)s^ 



but 

A* PA = N = Z " 1 
s^ s^ 

tho varianco-covarianco matrix. 

Thus 



(36) n - u R H = (Z - CX*)' (CSC)* 1 (Z - CX*) 
b R 0 b 



which is distributed as n-upC • The hypothesis can be rojoctod if 
the computed value of this quantity exceeds the tabulated value of 



**b, n-u,o< • 
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If a theoretical model has parameters X A the hypothesis 



that is 



H q : X - X‘ 



X 1 “ X 1 



x 2 - x 2 



x i “ x i 

The matrices C and W are 

u c n “ 1 and u^l " * T 
thus the statistic to be tested is 

(X* - X T )« Qp* 1 (X* - X T ) - % 
us^ u 

This is a more general derivation of the statistic offered earlier. 

To test a hypothesis concerning the value of a single parameter, 
regardless of the values which the other parameters have, the same pro- 
cedure would be followed. 



: x, - x 



1 



C - (1, 0, 0 



0 ) 



W - x. 



Then if the hypothesis is true 



T\ 2 



(n - u) Rjj - (xj* - x^ 1 ) 

-2 



1 



where Sj is the variance of Xj, the statistic (n - u) 






is distributed as Fj, This is simply the square of the 
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Student's t for (n - u) degrees of freedom. 

For a more detailed discussion o£ the matrix approach the reader 
is referred to chapter 4 of Hamilton (1964), 

To illustrate the matrix approach, the data of Hatch (1964) 
will be used. 



N - 



15.00 


0.00 


-0.10 


-0.01 


0.01 


-0.20 


-0.03 


0.00 


7.66 


0.13 


-0.20 


-0.01 


0.02 


0.03 


-0.10 


0.13 


7.42 


0.05 


0.00 


0.02 


0.01 


-0.01 


-0.20 


0.05 


7.31 


0.00 


0.03 


0.00 


0.01 


-0.01 


0.00 


0.00 


0.03 


0.00 


0.00 


-0.20 


0.02 


0.02 


0.03 


0.00 


2.32 


0.00 


-0.03 


0.03 


0.01 


0.00 


0.00 


0.00 


0.27 



Other data are 



s - 2,44 
s 2 - 5.95 
u - 7 
n - 15 



X* - 



4.05 
-2.46 
5.71 
1.74 
1.00 
-5.30 
- 0.12 

For the sake of an example, assume that some theoretical consid- 
erations predict to be: 
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x T = 



3.66 

o.oo 

6.83 

0.00 

0.00 

0.36 

0.00 



We wish to test the hypothesis 



The computed values are 



(X* - X A )= 



H q : X = X T . 



0.39 

-2.46 

- 1.12 

1.74 

1.00 

- 5.66 

- 0.12 



S H = (X* - X T )' 5-1 (X* - X T ) = 158.79 = 3.81 
U US ^ 41.65 

From the tables for F (Appendix IV) 

F 7, 8, 5# = 3.50 

Since ^ F 7 8 5 $* the hypothesis may he rejected at 

u 

the 5 # significance level. Rejection at this level of significance is 
termed significant . Can the hypothesis also bo rojoctod at the Vjo 
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level? 



f7, 8 , 17. " 6 * 84 

S H ^ F 7 so the hypothesis may not be rejected at the 17. 



significance level. 



3.4 Complex Models 

In many cases a simple two dimensional linear mathematical 
structure will not properly represent the given data. It may be that 
a polynomial of increasing order 

Y -A + A + /W + A x i n 

will better represent the true structure. The matrix least squares 
solution in procedure is identical to the previous case. Using the 
solution matrix (b Q> b^ ... b a ) and the weight coefficient matrix (Q) 
statistical hypothesis may be tested using the t tests previously shown. 
The hypothesis 



H o : A m m A'* @2 m A' ***/A 

is tested by computing 




s oVqu s i 



and tested against 



fc (l -c*/ 2 )(n 

The hypothesis may also be tested by 



- u) 



F- b/ 



s <Aii 



These F tests serve to assess the significance of the additional 
reduction in the residual sum of squares achieved by fitting b's in the 
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particular order adopted. Tho order of fit is important. In a 
polynomial this ordor is fixod. In ether structures tho order is deter- 
mined by the oquations of tho mathematical structure. 

Then fitting polynomials Ostle (1963) recommends, 

"Rather than sook a bottor fit in terms of a 
higher degree polynomial (i.c. a degree greator 
than 2) it is probably better to cast about for 
some othor functional form to roprosont the data." 

An excellent example of regression analysis applied to a gravity 
problem is given in the article by H. Wolf (1965). In this paper Dr. 
Wolf uses hypothesis tests to dotormino tho systematic trend of gravity 
differences along tho Europoan Calibration lino. 

It is evident that the matrix approach discussed in Section 
3.31 can be easily applied to complex models. 

Application of hypothesis tosts to a non-linoar model raises 
another interesting problem. It cannot bo shown that tho non-linear 
least-squares solution will always converge to oven a local minimum 
value of the woightod sum of tho residuals, V' PV, This problem is dis- 
cussed by Hamilton (1964). Ho concludos, 

"If the estimated errors are small onough that 
tho functions aro truly linear over tho rango 
of several standard deviations in oach parameter, 
the methods of testing linear hypothesis . . , 
can be applied in tho same way." 

Whon individual paramotors are tosted, either by the t or F 
tests previously demonstrated, any co-varianco between tho parameters 
is not considered. If the parameters woro truly indopondont tho choice 
botween individual tosts, i.o., the t tost, and the simultaneous 
matrix tests would be ono of individual preference. 

If the co-varianco botween paramotors is largo, the simul- 
taneous matrix test should bo used. 
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Example: In his work with tho MRA-1 Tolluromoter, Hatch (196*0 

investigated tho mathematical structure 

T-B = bj + b 2 sin (2)/- (A + b^)) + bj sin 

100 

-t ^ Z(,yr ^ A + b 6^ + b 4 Bin <i>T ( A + b?) ). 

100 Too 



Aftor least squaros ad justmont, results for January 19 were: 



b(l)= 


2.78 


+ 


1.29 


b(2)= 


0.32 




1.83 


b(3)= 


10.09 


± 


1.82 


b(4)= 


-2.57 


± 


1.80 


b(5)= 


-0.31 


± 


89.63 


b(6)= 


23.02 


± 


1.45 


b(7)= 


-0.06 


± 


2.92 



Ho made 16 observations thus 

(n - u) = 16 - 7 =9 



degrees of freedom. 
To test 



is computed from 



H » > 4 - "° 

\ “ b i /s bi 

t* = 2.78 = 2.15 
1.29 

t 2 = 0^2 = 0.175 
1.83 

t, = 16.09 = 5.54 
J 1.82 




1x31 = 



1.80 



1.43 



from the tablo 
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t 5 = 0.31 = 0.0035 
t 6 = 23.02 = 15.90 
t 7 = 0.060= 0.0206 

* 9 , = ? ' ZG ' 

thus tho only parameters which arc si.grri.fi cant at tho 5^ lovol are 
and b£. 

At the 1$ level 

* 9 . 1 *" 3.25 

therefore b^, b^, and b^ are significant at this lovol. 

Those tests at a significance level of 1/5 givo statistical sup- 
port to the conclusions drawn by Hatch (1964) that the error can be rep- 
resented by 

E = bj + bo sin 4 tt(A + b^). 

J 100 

To go a step farther with tho data given by Hatch, the hypoth- 
esis 

H o : % = 0, i = 1 to 7 

was tested on each sot of observations using the t test. Table 3 shows 
the results (Hatch, 1964). 
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Table 3 



Test of H Q : = 0 



bate/ 

Parameter 


19 

Jan. 


"21-22 

Jan. 


8 — 

Feb. 


25 

March 


T5 

April 


n — 

April 


Kujnbor of 
Rejections 


B i 


A* 


A 


R 


R 


A 


R 


3 


B 2 


A 


A 


R 


A 


R 


R 


3 


B 3 


R 


A 


R 


R 


A 


R 


4 


b 4 


A 


A 


R 


A 


A 


A 


1 


B 5 


A 


A 


A 


A 


A 


A 


0 


B 6 


R 


A 


A 


A 


A 


R 


2 


B 7 


A 


A 


R 


A 


A 


A 


1 



* A signifios acceptance, R signifies rejection. 



If the parameters, which rejected the hypothesis more than once 

in the six sets are used, the formula for the cyclic zoro error would bo 

E = b, + b£ sin Zjr A bo sin 4 tt (A + b^) 

100 J 100 

As a comparison of tho individual t test and the simultaneous 
matrix test consider tho following data taken on 23 March by Hatch 
(1964). 



u = 7, 


s o = + 3.03. s Q 


B( 1 ) = 


3.66 + 


O .65 


B(2) = 


-1.67 + 


0.87 


B(3) = 


6.83 ± 


0.93 


B(4) = 


1.75 + 


0.88 


B(5) = 


21.15 + 


15.39 


B(6) = 


0.36 + 


1.32 


B(7) = 


1.63 ± 


4.79 
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25.00 


- 2.95 


- 0.26 


2.58 


-0.27 


- 0.05 


0.11 




-2.95 


13.97 


-2.82 


-2.23 


-0.02 


1.60 


0.17 




- 0.26 


-2.28 


11.62 


-1.58 


0.01 


0.18 


0.14 


N = 


2.58 


-2.23 


-1.58 


13.05 


0.01 


-1.92 


-0.11 




-0.27 


-0.02 


0.01 


0.01 


0.04 


-0.08 


-0.01 




-2.05 


1 .60 


0.18 


-1.92 


-0.08 


6.16 


-0.28 




0.11 


0.17 


0.14 


-0.14 


-0.01 


-0.28 


0.42 


The inverse 


matrix is 


given as: 












0.05 


0.01 


0.00 


-0.01 


0.32 


0.62 


-0.00 




0.01 


0,08 


0.02 


0.01 


0.05 


-0.02 


-0.00 




0.00 


0.02 


0.09 


0.01 


-0.01 


-0.01 


-0.04 


N _1 = 


-0.01 


0.01 


0.01 


0,03 


0.00 


0.02 


0.03 




0.32 


0.05 


-0.01 


0.00 


M 


0.45 


0.72 




0.02 


-0.02 


-0.01 


0.02 


0.45 


0.19 


0.14 




0.00 


-0.05 


-0.04 


0.03 


0.72 


0.14 


2,51 



From the previous t tests ono would expect that B(4), B(5)» 
and B(7) might be equal to zero, Constraining all values to bo equal 
to their test value, the C matrix is equal to I, and Z is equal to the 
test matrix, X^. The hypothesis 



where 



H : X= X 
o 

3.7 
-1.7 

6.8 
0.0 
0.0 

0.4 

0.0 
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will bo tostod. Using equation ( 36 ) 

% = (x*-x T Vu (x*-x T ) = 44.697 = 0.693 . 

u us (7)(9. 19) 

O 

Sine© S H ^ F 7 , 18 , 5 ^ = 2 * 58 . 

u 

the hypothesis is accoptablo. The values of B(4), B(5)* and B(7) are 
probably equal to zoro, which is the same conclusion reached by the t 
tost. The values given for B(l), B(2), B(3)» and B(6) are not neces- 
sarily the final values. The new structure, containing only those para- 
meters, must be solved by another least squaros adjustment. For the 
data used previously the results of this readjustment and the change 
from the values given by Hatch aro 





Valuo 


Chanft' 


b(i) 


3.45 


.21 


B(2) 


-1.80 


-.13 


B( 3 ) 


6.73 


.10 


B(6) 


0.28 


.16 



Statistical tests could also bo applied to these values. 

The previous examplo could also havo bocn tested using the "R 
factor" test discussod by Hamilton (1964). This is a tost of the ratio 

r q 

o • For the details the reader is roforrod to page 157 of Hamilton. 

R o 

To point up the difference botween the individual tests and a 
simultaneous test, the data of Laurila ( 1965 ), shown in Figure 6, Section 
3.3 of this thesis will bo used. 

For this example assume that some theoretical consideration 
predicts that the values are 
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C' = 0.206 
K» = 2.10 • 



Using tho t tost 



t K' = 2.10 = 2.06 

1.2 

t C» = 0.206 = 2.06 . 

0.10 



From tho tables, ^ 2 ^,^ 2.06; thoroforo, thoso valuos of K’ and C 1 

are acceptable. 

The values of K' and C' can then be used as tho tost values for 
a simultaneous F test. 



(X* - X T ) = 



O.454 

-4.70 



[ 0.454 - 4 . 7 ] 


" 25.0 


272.4” 


" 0.454 




_ 272.4 


3824.28 


-.-- 4 . 7 , - 



(2)(0.2129) 



F = 195 , 683 . 



Tho tabulated value of 18,5/6 2. 58, thus tho hypothesis can be re- 



jected although tho test valuos are acceptable on an individual basis* 
One would expect this, because the Q"*^ matrix indicates a strong corre- 
lation between parameters* 



0 
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CHAPTER 4 



LEVEL OF SIGNIFICANCE 

“Hie soloction of tho significance level to bo used in a 
statistical tost is a matter of judgoment and oxporionco on tho part of 
the geodosist. It may bo choson before tho data is analyzod. Basi- 
cally tho problem is; how largo should tho critical region for rojoction 
bo, or what is the risk of committing a l^pe I error which tho geodesist 
is willing to accept? Dixon (1957) states, 

"If it is a mattor of great concern when a truo 

hypothesis is rejected, ex should bo small 

if it is a matter of groat concern that a hypothe- 
sis bo rojectod if thoro is littlo evidence against 
it w© should uso a largo cx " . 

A convention followod by statisticians is tho following: If a 

hypothesis is rejected at cX = 5$ it is said to be significant . A 
hypothesis rejected at oC = Yp is said to bo highly significant . 

When chosing a significance lovol for a test, one must keep in 
mind that the acceptance of a hypothesis is favored ovor rojoction of an 
alternative by apy tost. If an cX. of 1$ is chosen, tho critical area 
for rejection is small, thus rejection is less likely than if a 5$ 
level of significance is choson. 

Tho choice between cX = 5$ and (X. - 1$ must bo dictated by the 
circumstances surrounding tho problem. 

An entiroly differont approach to hypothesis testing can be 
developed by using a slightly different procedure. In this method tho 
significance level is not pre-seloctod. The tost statistic is computed 
tjy one of the methods previously aiscussod. 

The probability (found in the tables) associated with tho valuo 
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yieldod by tho data is taken as an objective moasuro of tho dogroo of 
support that tho data londs to tho hypothosis. Taking this approach 
tho test of significance does not lead to tho accoptanco or rejoction of 
the hypothesos; it moroly moasuros the strongth of boliof. In many 
problems encountered in goodosy and photogrammotry this may bo the 
bettor approach. 

Examplo: As a final stop in his analysis Hatch (1964) rep- 

resented each of his regression coefficients by a linoar rogro3sion 

bj = o + Cj 

b 2 = l-U P + C 3 

b 6 = m 6 0 + c 6 

Whore e is the partial vapor pressuro, and P is tho bararaotric 
pressure. 

Efy- a least squares adjustment ho obtained tho following: 

(P. 66) 





Valuo 


Standard Error 


Ml 


41.53 


9.1 


C 1 


- 8.08 


2.0 


M 3 


-12.09 


40.4 


c 3 


359.69 


118.1 


m 6 


-144.82 


30.8 


c 6 


43.80 


6.8 



The t statistic can bo usod to determine tho significance of 



t - x 

H ' 



each co-efficient 
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t 


Dogroo of Freedom 


Degree of Support 


M 1 


4.57 


4 


98$ 


C 1 


4.03 


4 


98$ 


M 3 


.30 


4 


^ 50$ 


c 3 


3.04 


4 


95$ 


m 6 


4.71 


4 


99$ 


c 6 


6.45 


4 


99$ 



b 2 and 



It can be concluded that thoro is littlo relationship between 
P since is significant at a level of loss than 50$. 
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CHAPTER 5 
CONCLUSIONS 

This thesis has prosontod tho statistical theory basic to 
hypothesis tests. Some of tho statistics commonly used in hypothesis 
tests havo boen discussed, and thoir applications to tho problems of 
geodesy demonstrated. 

The tests most applicable to geodotic problems aro obviously 
tne ones which use a statistic that does not requiro the true variance 
in its computation. Urns the Student's t statistic and tho F statistic 
are more useful. 

The t statistic is usod to test tho moan of a population with 
an estimated standard error, s, and a finite number of observations, n, 
against some other mean value. 

t = (X -//„) n* 
s o 

Fisher's F statistic is given as 




It can be used to test hypothesis concerning the variance of two inde- 
pendent sets of data. 

From the adjustment procedures outlined by Dr. Uotila (1966), 
all tho values noodod to compute tests statistics aro available. In 
most cases computations of the tost statistic nood only bo carried to 
three significant figures. This precision can bo obtained on an 
ordinary slide rule. 

By using matrix mothods the dosirod test statistic can bo com- 
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puted, during the adjustment solution, by an cloctronic computer. Tho 
simultaneous matrix tost is preforrod over tho individual t tost in 
cases whoro thoro is a largo correlation botwcon paramotors. Tho porson 
analyzing tho data need only compare tho computed statistic with tho 
tabulated values to make the hypothesis tests. 

Tho purpose of using statistical hypothesis tests in goodosy and 
photogrammotry is to guido tho usor to tho bost conclusions based on 
the data analyzed. It should bo emphasized that failure to roject a 
hypothesis doos not moan tho hypothosis is truo. If, on tho basis of a 
tost, a hypothosis is rojectod tho statement can be mado that thero is 
evidonce, from the data analyzod, that tho hypothosis is not truo. 
Conclusions drawn with tho aid of statistical tests aro thus supported 
by probability theory, in addition to tho judgement and oxporience of 



the scientist 
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APPENDIX I 

CUMULATIVE NORMAL DISTRIBUTION 
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Appendix I 

Cumulative Normal Distribution 





Values of V = J 


■ V </>(.r) 'lx for .V = O.OOfO 0112.00 

- 30 




0.09 


A — > 
! 


0.00 


0.01 


0.02 


0.03 


0.01 


0.05 


0.00 


0.07 


0.08 


o.o 


.5000 


.50 10 


.5080 


.5120 


.5150 


.5100 


.5239 


.5270 


.5319 


.5350 


0.1 


.53 os 


.5 r,s 


.5178 


.5517 


.5557 


.5500 


.5030 


.5075 


.571 1 


.5753 


0.2 


.570:5 


.5s:;2 


.5871 


.5010 


.5018 


.5087 


.0020 


.0004 


.0103 


.01 11 


0.3 


.0170 


0217 


.0255 


.0203 


.0331 


.O.’MS 


.0100 


.0 1 13 


.0180 


.05 1 1 


0.4 


.055 1 


.0501 


.0028 


.000 1 


.0700 


.0730 


.0772 


.0808 


.0814 


.0870 


0.5 


.0015 


.0050 


.0085 


.7010 


,7951 


.7088 


,7123 


.7157 


.7190 


.7221 


0.0 


.7257 


.7201 


.7.T21 


.7357 


.7380 


.7122 


.7154 


.7180 


.7517 


.7519 


0.7 


.75 SO 


.70,11 


.7012 


.7073 


.7701 


.773 1 


.7701 


.7701 


.78*23 


.7852 


o.s 


,7<SS1 


.7010 


.7039 


.7007 


.7005 


.8023 


.8051 


.SO 78 


.8100 


.8133 


0.0 


.SI 50 


.S ISO 


.8212 


.8238 


.820 1 


.8280 


.8315 


.8340 


.8305 


.8389 


1.0 


.8113 


.8138 


.8101 


.8185 


.8508 


.8531 


.8554 


.8577 


.8500 


.8021 


1.1 


.so 1:5 


.8005 


.8080 


.8 70S 


.8720 


.8710 


.8770 


.8700 


.8810 


.8830 


1.2 


.8840 


.ssoo 


.SSS8 


.8007 


.8025 


.89 i 1 


.8002 


.8080 


.8007 


.0015 


1.3 


.00152 


.00 10 


.900,0 


.0082 


.0000 


.01 15 


.0131 


.0147 


.0102 


.0177 


1.4 


.0102 


.0207 


.9222 


.0230 


.0251 


.0205 


.0270 


.0292 


.0300 


.9319 


1.5 


.0:5:52 


.05 15 


.0.’557 


.0370 


.0382 


.0301 


.0400 


.9118 


.0120 


.0111 


1.0 


.0152 


,0!0:4 


.0171 


.0181 


.0105 


.0505 


.0515 


.9525 


.0535 


.0515 


1.7 


.9551 


.0501 


.057.5 


.0582 


.0501 


.0500 


.00)08 


.9010 


.0025 


.0033 


l.S 


.00 11 


.0010 


.00,50 


.000 1 


.0071 


.0078 


.9080 


.9003 


.0000 


.0700 


1.0 


.971:5 


.0710 


.0720 


.9732 


.0738 


.9711 


.9750 


.0750 


.0701 


.0707 


2.0 


.0772 


.0778 


.0783 


.9788 


.0703 


.0708 


.0803 


.0808 


.0812 


.0817 


2.1 


.OS 21 


.0820 


.0830 


.0831 


.0838 


.08 12 


,08 10 


.9850 


.085 l 


.0857 


2.2 


.0S01 


.0801 


.0808 


.0871 


.0875 


.0878 


.0881 


.988 1 


.OSS 7 


.0800 


2.3 


.oso:5 


.0800 


.9898 


.0001 


.000 1 


.0000 


.0000 


.9911 


.0013 


.0910 


2.1 


.09 IS 


.9920 


.9922 


.0025 


.0027 


.0920 


.0031 


.9932 


.093 1 


.9930 


2.5 


.00.58 


.0010 


.9011 


.0013 


.9015 


.0010 


.9018 


.9010 


.9951 


.0052 


2.0 


.0058 


.00.55 


.0050 


.01*57 


.0050 


.0000 


.0901 


.0002 


.0903 


.900 1 


2.7 


.0005 


.0000 


.0007 


.0008 


.0000 


.0070 


.9971 


.0972 


.9973 


.007 1 


2.8 


.0071 


.0075 


.0070 


.0077 


.0077 


.0078 


.0079 


.0979 


.0080 


.008 1 


2.0 


.0081 


.0982 


.0082 


.0083 


.008 1 


.008 1 


.9085 


.9985 


.0080 


.0080 



Taken from Appendix, Table I (Hamilton, 196^-) 
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APPENDIX II 

STUDENT’S t DISTRIBUTION 
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Usod in conjunction with problems in statistics, the table of 
the Studont's t distribution permits tho evaluation of deviations ox- 
prossod in terms of estimates of standard errors for samples of various 
sizes. A given estimate of standard error is divided into a difference 
or deviation, to obtain t as a basis for a test of significance. The 
table is entered with tho number of degrees of freedom determined for 
the problem. The tabular entry in the column is tho value of t 
associated with the probability level indicated at the top of the 
column. This level expresses the probability of obtaining a difference 
as large as the one obtained due to chance. (Arkin and Colton 1959). 
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Student's t Distribution 



Degrees of 
Freedom 


Probability 


0.50 


0.10 


0 05 


0 02 


0 01 


1 


1.000 


0.34 


12.71 


31.82 


63 66 


2 


0.810 


2 . 92 


4 30 


6 m 


9 92 


3 


,705 


2 35 


3.18 


4.54 


G.M 


4 


.741 


2.13 


2.78 


3.75 


4 60 


5 


.727 


2.02 


2.57 


3.36 


4.03 


G 


.718 


1 .94 


2 45 


3 14 


3 71 


7 


.711 


1 . 90 


2.36 


3.00 


3 50 


8 


.706 


1 .86 


2.31 


2.90 


3 36 


9 


.703 


1.83 


2.26 


2.82 


3 25 


10 


.700 


1.81 


2.23 


2.76 


3 17 


11 


.697 


1 .80 


2.20 


2.72 


3.11 


12 


. 695 


1.78 


2.18 


2 . 68 


3.00 


13 


.694 


1.77 


2.1G 


2 . 65 


3.01 


14 


.692 


1.76 


2.14 


2 . 62 


2.98 


15 


.691 


1.75 


2.13 


2 60 


2.95 


1G 


.690 


1 75 


2 12 


2.58 


2.92 


17 


.689 


1 .74 


2.11 


2.57 


2.90 


18 


.688 


1 73 


2.10 


2 . 55 


2.8S 


19 


.688 


1.73 


2.09 


2.54 


2.80 


20 


.687 


1.72 


2.09 


2.53 


2.84 


21 


.686 


1.72 


2.08 


2.52 


2.83 


22 


.686 


1 .72 


2 07 


2.51 


2.82 


23 


.685 


1 71 


2.07 


2.50 


2.81 


24 


. 685 


1.71 


2.06 


2.49 


2.80 


25 


.684 


1 71 


2.06 


2.48 


2.79 


26 


.684 


1 71 


2 06 


2.48 


2.78 


27 


.684 


1.70 


2.05 


2.47 


2.77 


28 


.683 


1 70 


2.05 


2.47 


2.76 


29 


.683 


1.70 


2 04 


2.46 


2.76 


30 


.GS3 


1.70 


2.04 


2.46 


2.75 


35 


.GS2 


1 69 


2 03 


2 44 


2.72 


40 


681 


1 . 68 


2.02 


2.42 


2 71 


45 


.680 


1 .68 


2.02 


2.41 


2 69 


50 


.679 


1 .68 


2.01 


2 . 40 


2 GS 


00 


.678 


1 67 


2.00 


2.39 


2.66 


70 


. 678 


1 .67 


2.00 


2 . 38 


2.65 


SO 


.677 


1 . 66 


1 .99 


2 38 


2.04 


90 


.677 


1 . 66 


1 . 99 


2.37 


2.63 


100 


.677 


1 . 66 


1.98 


2 36 


2 63 


125 


.676 


1.06 


1 98 


2 36 


2 62 


150 


. 676 


1.60 


1.98 


2 35 


2 61 


200 


.675 


1 65 


1.97 


2 35 


2 60 


300 


.675 


1 . 65 


1.97 


2 34 


2 59 


400 


.G75 


1 65 


1 97 


2 34 


2 59 


500 


.674 


1.65 


1 . 96 


2.33 


2 59 


1000 


.674 


1 .65 


1 .96 


2 33 


2.58 

) 


00 


.G74 


1 .04 


1 96 


2 33 


2.58 



• The greater portion of this table taken from fl. A Fiaher'a "Statistieal Methods for ! tricart n 
Workers,” with the pcrmisuion of tho author and his publishers, Oliver and Boyd, London 



fiouree: Ilenrodurcd by permission from C, II. Goulden, Mtthods of Slotuhcal .inaly tt (N< w 
York: John Wiley &. Sons, 1939). 



Taken from (Arkin and Colton 1959) 
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APPENDIX III 
TABLE OF CHI-SQUARE 



61 



Tho tabic of Chi-squaro is entorod with the dogreos 
of froodom appropriate to tho problem. Tho row for the specified 
degroos of froodom is followed across to tho columns corresponding to 
o(/2 and 1 - o(/Z whoro tho thoorotical values of Chi-squaro needed for 
tho t©3 1 aro found. 
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Percentage Points of the x 2 Distribution* ’f 
Value . 14 of x 2 u.r>* where <\ is the probability that \ ? exceeds \ 7 tl and 

j* <t'(\ 2 ) ( h 2 = 1 — nr t 



N. c* 
n \ 


0.995 


0.990 


0.975 


0.950 


0.500 


0.050 


0.025 


0.010 


0.005 


i 


0.00+ 0.00+ 0.00+ 0.00 + 


0.15 


3 81 


5.02 


0.03 


7.88 


2 


0.01 


0.02 


0.05 


0.K) 


1.39 


5.99 


7.38 


9.21 


10.00 


3 


0.07 


0.11 


0.22 


0.35 


2.37 


7.81 


9 35 


11.31 


1 2.8 1 


4 


0.21 


0.30 


0.18 


0.71 


3.3(5 


9.19 


11.1 1 


13.28 


1 1.8(5 


5 


0.11 


0.55 


0.83 


1.15 


4.35 


11.07 


12.83 


15.09 


10.75 


G 


0.G8 


0.87 


1.21 


1.04 


5.35 


12.59 


11 15 


10.81 


18.55 


7 


0.99 


1.21 


1 .09 


2.17 


(5.35 


14.07 


10.01 


18.48 


20.28 


8 


1.31 


1 .05 


2.18 


2.73 


7.3 1 


15.51 


17.53 


20.09 


21.90 


9 


1.73 


2.09 


2.70 


3.33 


8.34 


10.92 


19.02 


21.07 


23.59 


10 


2. 1 G 


2.50 


3.25 


3.91 


9.34 


18.31 


20.18 


23.21 


25.19 


11 


2. GO 


3.05 


3.82 


4.57 


10.31 


14.(58 


21.92 


21.72 


2(5.70 


12 


3.07 


3.57 


1.10 


5.23 


11.31 


21.03 


23.34 


20.22 


28.30 


13 


3.57 


1.11 


5.01 


5.89 


12.34 


22.3(5 


21.71 


27.09 


29.82 


1 1 


4.07 


1.00 


5.(53 


0.57 


13.31 


23.08 


20.12 


29.14 


31.32 


15 


4. GO 


5.23 


ti.27 


7.20 


14.3 1 


25.00 


27.49 


30.58 


32.89 


10 


5.14 


5.-81 


0.91 


7.90 


15.34 


2(5.30 


28.85 


32.00 


31.27 


17 


5.70 


0.41 


7.50 


8.(57 


10.31 


27.59 


30.19 


33.11 


35.72 


18 


G.2G 


7.01 


8.23 


9.39 


17.34 


28.87 


31.53 


34.81 


37.10 


19 


0.84 


7.03 


8.91 


10.12 


18.31 


30.14 


32.85 


30.19 


38.58 


20 


7.43 


8.20 


9.59 


10.85 


19.31 


31.41 


34.17 


37.57 


40.00 


25 


10.52 


11.52 


13.12 


14.(51 


24.31 


37.(55 


40.05 


44.31 


40.93 


30 


13.79 


14.95 


1(5.79 


18.49 


29.3 1 


43.77 


4(5.98 


50.89 


53.07 


•40 


20.71 


22.10 


21 13 


20.51 


39.3 1 


55.70 


59.3 1 


03.09 


00.77 


50 


27.99 


29.71 


32.30 


31.70 


49.33 


(57.50 


71.42 


70.15 


79. 19 


GO 


35.53 


37.18 


40.18 


13.19 


59.33 


79.08 


83.30 


88.38 


91.95 


70 


43.28 


15 11 


•18.70 


51.71 


09.33 


90.53 


05.02 


100.12 


104.22 


SO 


51.17 


53.5 1 


57.15 


00.39 


79.33 


101.88 


100.(53 


112.33 


11(5.32 


90 


59.20 


01 .75 


(55.05 


(59.13 


89.33 


113.14 


118.11 


121.12 


128.30 


100 


07.33 


70.0(5 


74.22 


77.93 


99.33 


124.34 


129.5(5 


135.81 


1 10.17 



Adapted from the tables prepared by Catherine M. Thompson for Itiometrika, 
V() b 32; r(‘produced with permission of the editors of JUomdiiLu. 

t hor more than 100 degrees of freedom, percentage points x 2 »a of the x 2 distribution 
h(* obtained from the two-tailed percentage points A> of the normal distribution 
hv (he approximate relation, x 7 n. a n - b (2a) l/2 A>, with <* = i\ 
f or is thus the probability in one tail of the distribution. 



Taken from Appendix Table III (Hamilton, 1964) 
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APPENDIX IV 
TABLE OF F 
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“In tho tablo of F tho distribution of F is tabulated for tho 5/^ 
and 1$ rosidual lovols, Tho 5$ level value of F is indicated in ordinary 
typo, and tho 1$ lovol figures aro printed in bold faco typo. Entry in- 
to tho tablo is accomplished by rofcronco to tho appropriate column for 
tho degroos of freedom associated with tho greater varianco, and to the 
appropriate row for tho dogroes of froedom associated with tho smallor 
varianco. If tho calculated ratio botwocn the two variances (F) excoods 
tho valuo for F indicated in tho body of tho tablo for the 5$ level, 
there aro fewer than 5 chancos in 100 that the disparity botweon the 
calculated variances is duo to chanco; if F exceeds that recorded for 
the 1$ lovel, the probability is less than 1 in 100 that tho difference 
is accidental.''' (Arkin and Colton, 1959). 
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Tmux A-7 c. I’kkiTmiuis ok the /•’(»•», rd I ) i.vnuut T1U\ with Dhohki 
uk I’iuxuum !•, run Tin; Numkiiatih* and »> 2 run Tin. Dj.numin vmu 
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